Look-ahead Levinson- and Schur-type Recurrences in the Padé Table∗
نویسندگان
چکیده
For computing Padé approximants, we present presumably stable recursive algorithms that follow two adjacent rows of the Padé table and generalize the well-known classical Levinson and Schur recurrences to the case of a nonnormal Padé table. Singular blocks in the table are crossed by look-ahead steps. Ill-conditioned Padé approximants are skipped also. If the size of these lookahead steps is bounded, the recursive computation of an (m,n) Padé approximant with either the look-ahead Levinson or the look-ahead Schur algorithm requires O(n2) operations. With recursive doubling and fast polynomial multiplication, the cost of the look-ahead Schur algorithm can be reduced to O(n log n).
منابع مشابه
A look-ahead Schur Algorithm
The classical Schur algorithm computes the LDL T factorization of a symmetric Toeplitz matrix in O(n 2) operations, but requires that all the principal minors of the matrix be nonsingular. Look-ahead schemes have been proposed to deal with matrices that have exactly singular principal minors 9], 11]. Unfortunately, these algorithms cannot be extended to matrices that have ill-conditioned princi...
متن کاملNORTH- HOLLAND High Performance Algorithms for Toeplitz and Block Toeplitz Matrices
In this paper, we present several high performance variants of the classical Schur algorithm to factor various Toeplitz matrices. For positive definite block Toeplitz matrices, we show how hyperbolic Householder transformations may be blocked to yield a block Schur algorithm. This algorithm uses BLAS3 primitives and makes efficient use of a memory hierarchy. We present three algorithms for inde...
متن کاملA Block Toeplitz Look - Ahead Schur
This paper gives a look-ahead Schur algorithm for nding the symmetric factorization of a Hermitian block Toeplitz matrix. The method is based on matrix operations and does not require any relations with orthogonal polynomials. The simplicity of the matrix based approach ought to shed new light on other issues such as parallelism and numerical stability.
متن کاملA Block Toeplitz Look-ahead Schur Algorithm
This paper gives a look-ahead Schur algorithm for nding the symmetric factorization of a Hermitian block Toeplitz matrix. The method is based on matrix operations and does not require any relations with orthogonal polynomials. The simplicity of the matrix based approach ought to shed new light on other issues such as parallelism and numerical stability.
متن کاملLook-Ahead Procedures for Lanczos-Type Product Methods Based on Three-Term Lanczos Recurrences
Lanczos-type product methods for solving large sparse non-Hermitian linear systems have as residual polynomials either the squares of the Lanczos polynomials or the products of the latter with another sequence of polynomials, which is normally chosen to enforce some local minimization of the residual norm. In either case, these methods inherit from the underlying Lanczos process the danger of b...
متن کامل